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ABSTRACT 


The induced environment is the main source of contamination on a 
spacecraft. The contamination of surfaces by Volatile Condensable Ma- 
terial (VCM) present on spacecraft is being studied by many investigators 
both in the laboratory and space. In addition to VCM, the space environ- 
ment can also be a contributing factor in the contamination of surfaces. 
We have been investigating the effects of the interaction of the solar 
plasma with surfaces. Plasma particles impact at energies in the multi- 
kilovolt range and penetrate and become trapped in surfaces. Results 
show that protons can penetrate up to 600 monolayers into surfaces and 
reach saturation concentrations of 6%. 

1.0 INTRODUCTION 

Many investigators are conducting fundamental studies both in the 
laboratory! >2*’ anc j spaced'® on VCM contamination of surfaces as a 
function of temperature. Measurements are being made with the Tempera- 
ture-Controlled Quartz Crystal Microbalance (TQCM) and include quali- 
tative analysis of contaminant species, adsorption and desorption rates, 
sticking coefficients and changes in optical properties of surfaces. 

The TQCM? was developed in 1973 by Faraday Laboratories for NASA to 
monitor contamination between -60 and +100°C in space chambers. Re- 
cently, we completed construction of an extended temperature range 
TQCM system for space flight having higher mass sensitivity and ground 
command of temperature. Operating temperatures are attained either with 
radiative coolers or thermoelectric devices. The Air Force Materials 
Laboratory has a radiatively-cooled TQCM in flight on SCATHA** that was 
launched on January 30, 1979. JPL has a thermoelect rically controlled 
system ready for flight on a NOAA operational environmental satellite^. 
MSFC has both radiatively and thermoelectrical ly controlled systems 
awaiting flight on Shuttle^. The SCATHA system is shown in Fig. 1. 

It consists of a controller and two radiatively cooled sensor heads. 

The sensors are high mass sensitive 15-Mlz optically-polished quartz 
crystals plated with aluminum. 

The Principal Investigator for the SCATHA experiment is D.F. Hall. 

He has reported important data** on the measurement of surface contamina- 
tion enhanced by spacecraft charging in the multikilovolt range!!. In the 

* Supported in part by ONR Contract N00014-78-C-0373. 

For presentation at the 11th Space Simulation Conference, Johnson 
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experiment, one TQCM views space. The other periodically views space « 

and the sun at 1 rpm. Data from the sun viewing TQCM is significantly j 

different from the TQCM pointed into space. The mass accumulation rate 
of the sun viewing TQCM is relatively insensitive to sensor temperatures 
between -34 and +58°C. The space viewing TQCM temperature showed mass 
accumulation only for temperatures below -30°C and mass desorption at 
higher temperatures. The most probable reason for the difference is 
solar UV radiation modifying the contamination and increasing its de- 
sorption energy . Another effect that is independent of temperature 
and results in mass accumulation is the penetration and trapping of solar 
plasma into the aluminum surface plated onto the TQCM. We have been in- 
vestigating in the laboratory the effects impacts of protons and alpha 
particles present in the solar plasma have on surfaces. 

2.0 LABORATORY MEASUREMENTS 

A photograph of the plasma generator used for the measurements is 
shown in Fig. 2. A plasma is generated by first creating a narrow-energy 
ior. beam. Gas is fed into a flat cylindrical chamber with circular aper- 
tures on the axis. Electrons from a filament are made to spiral along 
this axis by an axially directed magnetic field generated by a solenoid. 

The spiraling electrons collide with the neutral atoms to produce ioniza- 
tion. The ions are extracted from the chamber (anode), pass through the 
focusing system, a series of metal tubes separated with ceramic insula- 
tors, and accelerate towards the target (cathode). The target is an alu- 
minum plated 20 MHz quartz crystal for measuring mass gain or loss during 
bombardment. A 20 MHz crystal was used because of its high-mass sensi- 
tivity of 8.78 x lO'l^g/cm^Hz needed to facilitate the measurements. The 
ion impact energy is controlled by varying the target potential. 

Beam fluxes of about i0l4 ions per second per cm2 are easily ob- 
tained with any gas using this beam generator. In our experiments, only 
noble gas atoms were chosen as bombarding ions for the following reasons. 

Monatomic beams of singly charged ions of noble gases can readily be 
formed without mass analysis. They are available in pure quantities. 

The range of atomic number and weight of these gases is more than adequate 
for these investigations. They are inert so that reactions with surfaces 
are minimized. The generator is operated in a Vac Ion pumped chamber. 

The effects of bombarding A1 with 0.35 mA He + at 2.5 keV is shown 
in Fig. 3. During the first 35 min of impact the target increases in 
mass because the mass of He + trapped in the surface was greater than the 
mass of A1 eroded by sputtering. For the next 15 min, the mass gained 
and lost is about equal, and the net mass change is approximately zero. 

After 50 min, the target became saturated with He. The mass loss rate 
is then equal to the sputtering yield of aluminum in He + . 

3.0 ION- SURFACE INTERACTIONS 

3.1 Surface Mass Loading 

The penetration and trapping of ions in a solid surface shown in 
Fig. 3 can be described by the mass rate of change per unit area, 
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dM/dt, of a target under bombardment. 

dM/dt = am^i - 3(t)m^ - un^i (1) 

The absorption coefficient, a, is the fraction of incident particles 
that penetrate and become trapped in the target. The beam flux is i in 
ions/cm2-s. The mass of the bombarding ion and of a target atom are, 
respectively, mj and m 2 . The re-emission coefficient, 3(t) is the number 
of trapped particles re-emitted as succeeding layers of target atoms are 
sputtered. The sputtering yield, u, is the number of target atoms e- 
jected per incident ion. At t = 0, 3(t) = 0, and as t becomes large, 

3(t) -ai. Once 3(t) =ai, the target will lose mass at a constant rate 
due to sputtering. 

3.2 Ion Penetration into Surfaces 

To penetrate beyond the first layer of atoms of a solid the inequality 

an 2/3 * 1 (2) 

must hold, otherwise the particle will probably be reflected. The atomic 
number density of the solid is n, and n^/3 represents the area number 
density seen by the incident particle with a total binary cross section, 

0 . Tins inequality is satisfied for protons and alpha particles in the 
solar plasma. 

Although the cross sections at these low energies are not known ac- 
curate./, Bohr^ has shown a reasonable cross section to be 


2rTa o E R Z 1 Z 2 m l + m 2 1 

2.72 (z 2/3 + z 2/3-jl/2 m 2 


(3) 


which is replaced by 

0 = ttR 2 (4) 

2 2/3 2/3 

whenever 0 »ua /(Z. + 1J ). At the lowest energies the binary col- 

lision cross section approaches the gas kinetic value. R is the distance 
of closest approach obtained as the solution of 


R m 2 E / (mj + m 2 ) = 2Z 1 Z 2 E R a 0 


exp [ -R( Z f 3 . /.„] 
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Subscripts 1 and 2 designate incident and target particles, respectively, 
with E the kinetic energy, m the mass, Z the atomic nunber, ao the Bohr 
radius and Er the Tydberg energy. 

To obtain the distribution of the penetrating ions it is necessary 
to consider their motion in the solid. Although the ions are probably / 1 

neutralized upon entering the solid this is expected to be of no con- 
sequency since hard sphere type collisions are assumed to occur. For 
protons and alpha particles, m^/m^ > 1, large angle scattering occurs 
at each collision giving an erratic rectilinear path which resembles a 
diffusion process. In fact, the slowing down process is treated like 
that of neutrons where the particle slowing down density for a plane 
flux i is 


q = 




( 6 ) 


where R is the "Fermi age." The density q represents the number of par- 
ticles per cm3 per second that slow down past an energy F at a given 
depth x in the solid. If n^/m^ > 1 the angular deflection of the inci- 
dent particle from its direct ; on of travel is small and the particle 
straggles about the penetration depth with an approximate Gaussian dis- 
tribution. 

Because the notion of range is complex, it is customary to define 
a linear, a vector and a projected range. The linear vvnge is the total 
rectilinear path traversed by the particle to the point of stopping while 
the vector range is the vector distance from point of entry into the 
solid to the end point. Projection of the vector distance along a 'ine 
normal to the surface gives the projected range. An approximate relation 
between the average linear range [Ri ] and the average projected range 
[Rp] (or mean range) was developed by Linhard and Scharffl3 as 


1 r l 1 “ tty <* * 1 sf> 


( 7 ) 


for an inverse square potential. For n^/mj ^ 1, approximation (7) is 
quite accurate but, for m2/m^ > 1 the approximation has application to 
a portion of the low energy region only. 

Scattering near the end point of the particle travel is probably 
a mar.y-body collision problem and not the binary one assumed to obtain 
the distributions discussed in the previous section. If relation (2) 
is used as a necessary criterion for halting the motion of the incident 
particle, the cross section must be about the gas kinetic value, 
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a « 5 x KT^cm^; cross sections of these magnitudes correspond to en- 
ergies the de Brolie wave length is comparable to the lattice separation 
indicating that more than one target atom contributes to a collision. 

Andersen and Sigmund 14 have presented a model for non-binary scat- 
tering which is of importance at much higher energies than indicated by 
(2). It recognized the fact that, with such small mean free paths, cor- 
relations between successive collisions cannot be neglected in the low 
energy region. The model assumed involves a set of crystal lattice 
rows which are resolved into stacked symmetric rings of atoms. The in- 
cident particle moves along the symmetry axis of these rings with the 
projectile-ring interaction as the basic event with the ring presenting 
an effective potential field to the incident atom. 

If the incident particle is aligned along the symmetry axis, close 
collisions become improbable because the field keeps the particle near 
the symmetry axis to produce channeling. For the interpretation of pene- 
tration distributions known in the literature as "supertails" it is 
important to investigate the low energy collisions to determine the den- 
sity of vacancies, interstitials and self- interstitials created. From 
this, estimates of further projectile penetration by diffusion after 
the slowing down process can be made. Using a Bom-Mayer potential for 
each lattice atom, a new reflecting (or stopping) energy replaces the 
one determined through (2) and is found to be 

R refl * P Ae ' L/a W 

whenever mj » m?. Constants p and L define the number of particles in 
the ring and radious of the ring element; A and a are the potential 
parameters. Reflection energies given by (8) are somewhat higher than 
that found through (2); for alpha particles in heavy materials, typical 
estimates are 5 10 eV. 

3.3 Diffusion of Trapped Particles 

Channeling explains the penetration to ranges of a few microns for 
particles in the keV range but "supertails" extend to depths of several 
microns. A possible explanation to these deep penetrations is the steady 
state diffusion approach of Sparks • He assumes the diffusion of vacan- 
cies and self-interstitials created by the incident particles is unaf- 
fected by the relatively small number of implanted ions. With the bom- 
barding particles forming substitutionals through the annihilation of 
vacancies two general forms of the integral penetration distributions, 
%(x), result which may be expressed as 

n s (x) » expf-k(x-l)] ; x > l (9) 
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n (x) * - ; x > l . (10) 

5 (x - l + d) 2 r 

The constant d depends upon the vacancy density interstitial density 
and the diffusion coefficient; l is located approximately at the edge 
of the stopping region. Distribution (10) normally occurs if vacancies 
and self- interstitials are created and annihilated in pairs. For dis- 
tribution (9) to obtain the vacancies, the self- interstitials must be 
created at different rates which means that most of the interstitials 
must be of the bombarding particle type. 

4.0 ION-BOMBARDMENT OF ALUMINUM 

Laboratory measurements were made on aluminum to determin. the 
effect the ions present in the solar plasma have on low atomic weight 
surfaces. The main constituents of the plasma are proton and alpha 
particleslb having energies in the multilcilovolt range. Heavier ions, 
such as iron, were not considered because their densities are too low 
to produce measurable effects. 

Once the protons penetrate the surface, little will diffuse out 
because of the affinity of aluminum for hydrogen^-". Any alpha particles 
that penetrate the surface will readily diffuse out because of the inert 
nature of helium. The principle effect of alpha impacts will be that 
of sputtering the surface. 

To determine the saturation concentration of protons for aluminum, 
the distribution function q given by (6) can be used. The distribution 
is for the case mj > m 2 - It has a maximum near the surface because of 
large angle scattering that occurs for the case in which the incident 
particle mass mj is less than the target atom mass m 2 * 

The mean penetration depth [FL] of the particles needed to scale 
the distribution function can be found from the p - e relationship of 
Lindhard and Scharff,18 which relates the range-energy of particles 
in solids: 


„ - d n A 2 -i 166 A 1 

0 V1 ^777 


e = 


33 


Z 1 Z 2 (Z 1 /3 + z 2 /3 ) 1/2 A x ♦ A 2 


( 11 ) 


( 12 ) 


193 



where the subscripts 1 and 2 refer to the incident ion and the target 
atom, and is the energy of impact in keV. 

Davie's group has shown that the p - e relationship is accurate to 
about 20$ for e < 10. Values for p, knowing e, have been plotted by 
them for 0 < e < 10 and 0 < p < 100.151 

The absorption coefficient a was found, experimentally, to be about 

0 . 01 . 

The velocity v at which the surface is sputtered was calculated 
using a yield, for A1 in He + is the same as that for He ++ (alpha par- 
ticles). Experimentally, the yield for He + at 8 keV was found to be 
0,3. The average erosion velocity of A1 using these yields for He + and 
He ++ is v» 10'°cm/hr. By assuming that diffusion is zero, the satura- 
tion concentration of trapped hydrogen in aluminum was found. 

The saturation distribution of hydrogen in aluminum, in hydrogen 
atoms per aluminum atom N^/Nyy , is shown in Fig. 4. 

5.0 CONCLUSIONS 

A spacecraft traveling beyond the magnetosphere passes into the 
tenuous upper atmosphere of the sun known as the solar plasma. The 
plasma is composed mainly of electrons, protons, and alpha particles. 
Protons and alpha particles will impact at energies ranging up to 10 
keV during periods of a quiet run and approximately 20 keV during solar 
activity .lb 

Low atomic weight metal surfaces will be damaged by the penetration 
and trapping of protons. Surface properties would be affected by the 
formation of hydrides in concentrations approaching 61. It will take 
several years before surface properties would be significantly affected 
because of the low density of protons in the solar plasma. Damage 
would be accelerated during solar activity and for spacecraft passing 
near the sun. 

Effects on nonconducting materials such as plastics, which are 
composed of large amounts of hydrogen, will be more pronounced because 
most of the plasma will be absorbed. For these materials, changes 
would become apparent in a few years. The protons and alpha particles 
that penetrated these surfaces will displace surface atoms. It takes 
about 25 eV to displace an atom from its lattice pointed. In a year 
the plasma will produce about 10 17 atomic displacements. In surface 
coatings, the displacements will produce discoloration through the 
creation of color centers and effect the solar absorptance. 


194 


I 



6.0 REFERENCES 


1. Muscari, J.A., and Jacobs, S., "Silver-Teflon Contamination UV Radia- 
tion Study", Proceedings USAF/NASA International Spacecraft Contami- 
nation Conference, AFNC-Tk-78-190, NaSA-cI)- 2039 (1978) pp. 111Z - ~ 


2. Richmond, R.G. and Kelso, R.M., Effectiveness of Shuttle Orbiter 
Payload Bay Liner as a Barrier to Molecular Contamination from Hy- 
draulic Fluids", Proceedings USAF/NASA International Spacecraft Con- 
tamination Conference, AFML-Tft-78-190, NASA-CD-2039 (1978) pp. 846 - 
WT. 

3. Visentine, J.T., Richmond, R.G., and Kelso, R.M. , "Experiment to 
Measure Molecular Outgassing Rates from Shuttle Orbiter Flexible 
Reusable Surface Insulation (FRSI)", Progress in Astronautics and 
Aeronautics, Vol. 56, AIAA (1977) pp. 197 - 213. 

4. Hall, D.F., "Flight Experiment to Measure Contamination Enhancement 
by Spacecraft Charging", Proceedings SPIE Conference, North Holly- 
wood, CA (1980) Vol. 216 

5. Miller, F.R., "Shuttle Induced Environment Contamination Monitor", 
Proceedings USAF/NASA International Spacecraft Contamination Con- 
ference, AFML-TR-78-190, NASA-CD-2039 (1978) pp. 534 - 566. 

6. Lynch, J.T., "Quartz Crystal Microbalance (QCM) Monitor of Contami- 

nation for LF.S-8/9", 9th Space Simulation Conference, NASA CP-2007 
(1977) pp. 211 - 219. " 

7. McKeown, D., Corbin, W,E. Jr. and Naumann, R.J., "Thermoelectrically 
Cooled Quartz Crystal Microbalance", 7th Conference on Space Simu- 
lation, NASA SP-336 (1973) pp. 345 - 35F! 

8. McKeown, D., Claysmith, C.R. , Breckenridge, W.T. Jr., Duraner, R.S., 
"Temperature Controlled Quartz Crystal Microbalance System", Report 
AFML-TR-77-83, AFML/MBE , Wright -Patterson AFB, OH (1977). 

9. Maag, C.R.,"N0AA-B Contamination Monitoring Instrumentation", Pro- 
ceedings SPIE Conference, North Hollywood, CA (1980) Vol. 216 - 1$. 

10. McKeown, D. and Claysmith, C.R. , "Quartz Crystal Microbalance Systems 
for Shuttle Contamination Measurements", Proceedings USAF/NASA Inter- 
national Spacecraft Contamination Conference, AFML-TR-78-lyO, NAsA- 
Cfrim (1*>?8) pp. 60S - 629. 


195 


11. Kuch, G.A. and Holaday, B.H., "A Description of the Space Test Pro- 
gram P78-2 Spacecraft and Payloads" Report SAMSO TR-78-24 (1978). 

12. Bohr, N., Penetration of Atomic Particles Through Matter , Det. Kgl. 
Dan. Vidensk, Selsk., Mat. - Frs. Medd., (1948)7 

13. Linhard, J., Scharff, M. , and Schiott, H.E., Range Concepts and 
Heavy Ion Ranges, Dt. Kgl. Dan. Vidensk. Selsk. , Mat. - Fys. Medd., 
Vol. 33, No. 14 (1963), pp. 643 - 721. 

14. Andersen, H.H., and Sigmund, P., "Simple Non-Binary Scattering 
Model Applicable to Atomic Collisions and Crystals at Low Energies", 
Dt. Kgl. Dan. Vidensk. Selsk., Mat. - Fys. Medd., Vol. 4 (1966), 
pp. 723 - 752. 

15. Sparks, M. , "Theory of Supertails of Ions Bombarded into Crystals", 
Phys. Rev. Letters, Vol. 17. (1966), pp. 1247 - 1249. 

16. Freeman, J.W. , Van Allen, .T.A., and Cahil, L.J., "Explorer 12 Ob- 
servations of the Magnetospheric Boundary and the Associated Solar 
Plasma on Sept. 13, 1961", Beophys. Res., 68, (1963), pp. 2121 - 
2130. 

17. Bergstrom, I., Brown, F., Davies, J.A., Geiger, J.S., Graham, R.L., 
and Kelly, R. , "On the Electromagnetic Separation Method of Preparing 
Radioactive Sources for Precision 3 - Spectroscopy", Nucl. Instr. 
Methods, 21, 249 (1963), pp. 249 - 261. 

18. Lindhard, J., and Scharff, "Energy Dissipation by Tons in the keV 
Region", Phys. Rev., 124 (1961), pp. 124 - 132. 

19. Davies, J.A., McIntyre, J.D., and Sims, G., "The Range of C^7 Ions 
of keV Energies in Germanium", Can. Jour. Chem., 40, (1962), pp. 

1605 - 1621. 

20. Seitz, F., and Koehler, J.S., "Displacement of Atoms During Irradi- 
ation", Solid State Physics, (Academic Press, New York, 1956), pp. 

307 - 44$; 


196 


1 




ORIGINAL PAGE IS 
OF POOR QUALITY 


« 


Fig. 1 SCATHA TQCM System without RPA Installed 
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Fig. 4 Saturation concentration of hydrogen m aluminum 
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